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Abstract 

We have investigated the predictions of neutrino oscillations within extended Anti-GUT, 
based on a large gauge group - the Standard Model and an B — L abelian group assumed 
with separate gauge fields coupling to each family of quarks and leptons - i= * 2 3 (SMGi x 
U(1) B -L,i)- We take into account corrections concerning a crude way of accounting for the 
number of ways the vacuum expectation values of the Higgs fields of the model can be ordered 
in the Feynman diagrams yielding the mass matrix elements. Performing these corrections 
in a balanced way between the charged fermion sector and the neutrino sector (case I) leads 
to the previous version of our model that was a fit to the MSW small angle solution region. 
The fit considered here is marginally worse that the previous fit inasmuch as it predicts a 
somewhat lower solar neutrino mass square difference (which might though be the only way 
of avoiding the day-night exclusion region). The other two versions (case II and case III) 
are characterised by the heaviest of the see-saw neutrino matrix elements that are composed 
either of v Rll and v Rt or two v Rt , These cases II and III fit the small mixing angle MSW 
scenario very well. 
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1 Introduction 



The latest results of Super-Kamiokande collaboration Q suggested that the day-night effect 
spectra disfavour the MSW (2| small mixing angle solution (MSW-SMA) at the level of 95% C.L. 
However, it might be that the MSW-SMA could not be excluded by the experiments because 
the measurement of the day-night effect is rather difficult. Therefore, we prefer to allow the 
possibility of the MSW-SMA, since we have recently predicted the MSW-SMA within extended 
Anti-GUT Q framework. 

It is also important to study whether the Anti-GUT model could, after all, predict the 
MSW large mixing solution (MSW-LMA): in our calculations we have used a technical correc- 
tion [|| ( "factorial factor correction" ) which takes into account the number of Feynman diagrams 
contributing to a given mass matrix element in a crude statistical way. In fact we used the pa- 
rameters from a fit to quarks and charged leptons mass matrices without "factorial factors 
correction", while we partly used them in the neutrino oscillation calculations. The major aim 
of the present article is to treat the charged fermions and the neutrinos on an equal footing with 
respect to this technical correction. We have assumed that one of the parameters is fixed to be 
unity as in earlier works |]]. In this article we will take into account that this parameter is not 
a priori unity. However, even these smaller variations do not bring our model to fit MSW-LMA 
but still only the MSW-SMA domain. 

Although our extended Anti-GUT model is rather restricted with respect to the choice of 
charge combinations for the various fermions and even for the Higgs fields introduced to break the 
gauge group assumed in order to keep the already achieved good fits, there are a few possibilities 
for playing around. One of the most important possibilities is that we can vary the quantum 
numbers for the Higgs field, called (j) B -L, the vacuum expectation value (VEV) of which is giving 
the scale of the see-saw particles (the right-handed neutrinos) by breaking the (B — L) charge 
(which is assumed to be gauged even in the form of one (B — L)i for each family where i denotes 
a generation) and thereby the overall scale of the neutrino oscillations. 

In the following section, we describe briefly the extended Anti-GUT model for the calculation 
of neutrino mass squared differences and the mixing angles. Then, in the next section, we put 
forward the mass matrices. In section 4 we describe the "factorial correction". Then in section 
5 we will present the experimental data and the "old" fits on the charged lepton and the quark, 
while we in section 6 give the results of our calculations. Section 7 contains our conclusion and 
resume. 

2 The extended Anti-GUT model 

The Anti-GUT model is characterised by a gauge group which for each family has a set of (family 
specific) gauge fields. In addition there is an extra U{1) gauge group called U(l)f. The details 
of the couplings of the latter is largely specified by the requirements of the model being free 
of gauge and mixed anomalies. So each generation gets its own system of gauge particles, i.e. 
gauge fields come in generations just like the fermions. The extension consists in providing each 
family also with a gauge field coupling to the (B — L)-charge of that family alone. That is to say 
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10 2 GeV » 10 12 GeV ^5 ~ 10 19 GeV 



Fig. 1 

The energy scale of the 
vacuum expectation val- 
GeV ues of the Higgs fields 

within the extended Anti- 
GUT model. 



we postulate for each family^ the subgroup of the well-known grand unification group, £0(10), 
which consists only of those generators that do not mix the different irreducible representation 
of the Standard Model, SMG, x U{l) B _ L>i = 517(3) x SU{2) x 17(1) x 17(1), where SMG denotes 
the Standard Model gauge group. The £7(1)/ that seems a strange extra ingredient in the "old" 
Anti-GUT model - although called for in order to be able to fit the quark masses || - is in the 
extended version just a linear combination of {B — L)i and yi/2 which happen to have separate 
anomaly cancellations for the right-handed neutrinos and for the rest (true Standard Model). 
Since the "old" £7(1)/ quantum charge is a linear combination of (B — 7)2, {B — 7)3, 2/2/2 and 
?/3/2 any anomaly constraint for £7(1)/ and the extended model charges must be automatically 
satisfied. Thus there formally can be no hindrance for having U (1) / as well as separate 77(l) B _ i 2 
and £7(l) fl _ £>3 gauge groups. However, such a possibility would imply that there would be a 
linear combination of the charges that would decouple from all the fermions and thus should 
not really be considered. In fact, it is a part of the arguments for the Anti-GUT model that 
one decides to ignore gauge fields which decouple from all observed fermions on the ground that 
they would have very little phenomenological relevance. 

The Higgs fields in our model are assigned quantum numbers under the gauge group "ex- 
tended Anti-GUT" = i= * 2 3 (SMGi x U{\) B _ L>i ) which are determined by seeking a fit to the 
quark and lepton spectra including its mixing angles. It should, however, be kept in mind that 
the possibilities for fitting these charge quantum numbers are few, since they are only allowed 
to take quantised values - analogous to that these quantum numbers for the quarks and leptons 
only take simple quantised values - and thus the model remains very predictive. The model ends 
up having only 7 Higgs fields falling into four classes according to the order of magnitude of the 
expectation values^: 

1 ) The smallest VEV Higgs field in this model plays the role of the Standard Model Weinberg- 
Salam Higgs field, 4>wsi 

with the VEV at the weak scale being 246 GeV/V2. 

2) The next smallest VEV Higgs field is also alone in its class and breaks the common B — L 
gauge group £7(1), common to the all the families. This symmetry is supposed to be broken 
(Higgsed) at the see-saw scale as needed for fitting the over all neutrino oscillation scale. 
This VEV is of the order of 10 12 GeV and called (f> B _ L . This field was not assumed to exist 

1 Really, we should say "proto-family" because (1) these proto-family get mixed before they are identified with 
the true families (2) even the right-handed charm- and top-quarks get permuted. 
2 The quantum numbers of the seven boson fields are shown in Table 1. 
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in the "old" Anti-GUT modelfl. 



3) The next 4 Higgs fields are called £, T, W, and x and have VEVs of the order of a factor 
10 to 50 under the Planck unit. That means that if intermediate propagators have scales 
given by the Planck scale, as we assume, they will give rise to suppression factors of the 
order 1/10 each time they are needed to cause a transition. The field x was introduced 
for the purpose of the study of neutrinos and was not present in the "old" model. 

4) The last one, with VEV of the order of the Planck scale, is the Higgs field S, which 
gives little or no suppression when it is applied. Therefore a transition amplitude in first 
approximation is not noticeably suppressed by this field S. Thus it gives rise to ambiguities 
in the model and its presence are not easily distinguished in phenomenology. Only if we 
take the VEV discernibly different from zero and/or make use of the "factorial corrections" , 
there is a possibility to observe phenomenological consequences of the field S. 

Therefore it is part of the model also that all physics beyond Standard Model is in un- 
extended Anti-GUT first a couple of orders of magnitude under the Planck scale and in the 
extended one begins at the see-saw scale. So there is pure Standard Model to see-saw scale 
and consequences of the Standard Model such as the GIM mechanism of no flavour changing 
neutral current are valid with the correspondingly very high accuracy, by far unaccessible to 
present days experiments. 

Since we have non-zero VEVs of scalar fields which are doublets under two different "proto"- 
generation specific SU(2)'s such as S, W, T, £ (while <fi ws has to be doublet for an odd number 
of generations), the mass eigenstates of the quarks and leptons do not belong to definite repre- 
sentations of the "proto" -generation specific 5'L r (2)'s. There is therefore no hindrance in having 
non-zero mixing angles, but they are of course suppressed corresponding to the need for the 
breaking VEVs being involved. 

In our previous article [H] we formulated the quantum numbers of our model with the use of 
the quantum number for the subgroup of the full gauge group called U(l)f, but it is more elegant 
without this abelian gauge group. The U(l)f quantum number, Qf, used in earlier articles - 
both in extended Anti-GUT (where it can be avoided), and in "old" Anti-GUT (where there are 
no right-handed neutrinos and thus no anomaly free (B — L)'s, so that U(l)f is unavoidable) - 
is related to the quantum numbers of Table 1 by 

Q / = (B-L) 3 - (jB -L) 2 + 2(|-|) . (1) 

The calculation of the mass matrices in our model consists in evaluating for each mass 
matrix element which of our seven Higgs fields are needed to provide the difference in quantum 
numbers between the right-handed Weyl components and the fermion in question to the left- 
handed ones. Then, one imagines that the propagator fermions all have Planck or fundamental 
mass scale masses in a Feynman diagram which is often a long chain of interactions with the 
successive Higgs fields that were needed. The order of magnitude of the diagram is then of the 
order of the Planck scale multiplied by a "suppression factor" for each Higgs field used. Now 
we assume that there are of order unity random couplings all along the chain and we therefor 
take for the value of the matrix element the product of the suppression factors times a random 

3 The influence, on the charged lepton mass matrix, of the new Higgs field, x> i s om y on the off-diagonal 
elements which are remain suppressed and do not dominate any masses or mixing angles. 
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Table 1: All U(l) quantum charges in extended Anti-GUT model. We presented here the three 
different possibilities of the (f) B - L quantum charges. The symbols for the fermions shall be 
considered to mean "proto" -particles. Non-abelian representations are given by a rule from the 
abelian ones, see after Eq. (2). 
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factor of order unity. These factors are taken as random numbers and at the end a logarithmic 
averaging of the resulting mass eigenvalues and mixing angles is taken. 

3 Mass matrices and Higgs quantum numbers in the extended 
Anti-GUT model 

To write down the mass matrices we need the quantum numbers of the Higgs fields. However, 
there is the freedom that we can modify the charge assignments without too much change in 
the predictions if we change the fields by adding to their quantum number assignment the 
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quantum numbers of S. Thus we shall take the point of view that we really do not know from 
the already developed phenomenology of the charged sector (quarks and charged leptons) the 
abelian quantum numbers except modulo those of S. 

To take into account the ambiguity of the choice of the quantum charges of Higgs fields, due 
to the VEV of S being of order of unity in Planck units, we could parametrise the quantum 
number combinations by integer parameters, a, /3, 7, 5, rj and e telling the number of additional 
S*-quantum number combinations counted from some starting value of the quantum number 
combinations of the Higgs fields which are given in Table This consideration leads to the 
following Higgs field quantum charges: 



W 



X 



case I 
case II 
case III 



1 2 2 , 
;,0,--,-,0) 



6'"' 3'3' 

6 3 2 3 3 ; 
1 12 2 

D O 3 3 

g,-g,0,-,--,0) +7 5 

0,0,0,0,-1,1) +rj s 
0,0,0,1,0,1) +eS 
0,0,0,0,1,1) +eS 
0,0,0,0,0,2) +eS . 



(2) 



Table 1 and Eq. (|2|) is simplified by only containing the abelian quantum numbers, but we 
do in fact imagine that in our model we have the proto-generation specified by the following 
rule: For each "proto" -generation i the non-abelian representations - of SU(2)i and SU(3)i - 
are the smallest dimension ones obeying the restriction 



ck/2 + + yi/2 = (mod 1) 



(3) 



Here is the triality being 1 for 3,-1 for 3 and for 1 or 8, respectively. The duality di is 
for integer spin SU(2)i representations while 1 for half-integer SU(2)i spin. 

Concerning the see-saw scale determining field (j) B -L, we have listed three "promising" choices 
of the quantum numbers in addition to the shifting S-quantum numbers choice. With these 
quantum number combination we get the parameterised mass matrices as follows, where we 
simply write S, W, etc. instead of (S), (W), etc. for the VEVs: 



the uct-quarks: 

U. 



V2 



the dsb-quarks: 



M D ~ 



V2 



v S 3 T+*(£t)3 

g-2-2a+40 — f-5 



g2+a-2l3-j + S^y\rp2^ g2a -/3- 7 + <5 (X^\ \ 2 
S 2+ a -2l3+S W \ T 2 S 2a-0+S ^ W U2 T 

gi+s Q,- 1+a+ (3 +S yy] yt 



(4) 



W 2 (T^C S- 



g-2-a+2f3—,-5-\y^jr\^2£ g2-3/3-~,-Srp3^ 
g-2- a +2l3-5]y^rpj^2 g2-30-Srp3 

WT 



3-2u+4f3-S 



W 2 (T^) 4 S 



(5) 



1-C-/3-S 



4 The starting values were though not the ones of Table 1 as seen from equation ([[2]) below. 
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the charged leptons: 

) 3 



M E - ■ ^ s 



V2 



(6) 

the Dirac neutrinos: 

(^l +a _2/5+2 7 +5^tj 1 2^t)2 g'2+ a -2i3+3j+S^y\rp2^-\^3 g2+c,-20+3-y+S-n (W^\T 2 U^^X^ 
g-2+ a +0— 1 +s+v\Y'tT^£x^ S~ 1+a+f3+s+ri W^T^x' ( S~ 1+a+l3+s W^T^ 

(7) 

Now we have to get the right-handed Majorana neutrino mass matrix to be able to calculate 
the effective neutrino mass matrix using the see-saw mechanism ||, ^, [TTJ. However, there are 
three different choices of the quantum charge of (ft B _ L , none of which give results that can be 
excluded immediately: the (1, 3)-component of the right-handed mass matrix is dominant (called 
case I), another choice is the (2, 3)-component is the dominant one (called II), and third one is 
the (3, 3)-component as the dominant one (called III). In our previous paper, the right-handed 
neutrino mass matrix is considered as case I. 

The right-handed Majorana neutrinos is given in the case I by: 

M R ^(<p B _ L )\ P^-'xW P+t-'tf (8) 



In the case II: 

In the case III: 
Mr 



/ S -2-2 J+v -c x ^2 S l +7+ „_e x t£ 51-HT-.A 

Mr ~ (</>b-l) I S^+^xH S—-x ] (9) 
51+7-^ s -* S~^x 



s-i-**>-<tf)H s^{x ] ? s^ x ] I • (10) 



Note that the quantum numbers of our different Higgs fields are not totally independent in 
so as far as there is a relation between the quantum numbers, 

Q x = 3Qw - 9Q r + (-6 - 3a + 9/9 + 77) Q s , (11) 

and thus the Higgs field combinations needed for a given transition are not unique, so the choice 
of the largest contribution for each matrix element must be selected. To compare with our earlier 
work (except for an d Table 1, we remark that the quantum number combination used in 
that work is obtained by putting 

a = (3 = 1 , 7 = 0, 6 = -1 , 77 = , e = 0. (12) 
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4 The number of orderings correction for the Feynman diagrams 



The technical detail called "factorial correction" consists in the following argument for a 
Feynman digram counting correction: the external lines signifying the attachment of the various 
Higgs fields used to make the quantum numbers match, can be put in several different sequences 
along the chain of fermion propagators in the diagram (see Fig. 2) providing the transition. 
For all these different sequences quite different fermion propagators are used. These differently 
ordered VEV-attachments are expected to be statistically independent and should be added with 
random phases. That means that one has a random walk in the complex plane (of amplitude 
values). The average of the numerical square of the amplitude for the mass matrix element, 
say, is getting additive contributions from the various diagrams and so goes up linearly with 
the number of diagrams. This means then that the amplitude goes as the square root of the 
number of diagrams with different orderings of the attachments of the different Higgs field 
symbols designating the action of these Higgs field VEV's. If we, for example, think of a matrix 
element for the electron mass, it turns out that in the model, it is necessary to use the vacuum 
expectation value for the Weinberg-Salam Higgs field {(f>ws) once, for the £ twice, for the T 
twice, for the W once, and for the S field n times where n depends on the specific assignment 
of the quantum numbers to the other fields. 

The number of orderings in which we can have the attachments of the VEV's for the electron 
mass generating diagram (in first approximation) is then (6 + n)!/(2! 2!n!). This means that we 
expect the amplitude statistically to be 4> ws ^ 2 T 2 WS n y/(6 + n)\/(2\ 2! n\). When these crude 
estimates of this correction were taken into account in the charged quark and lepton fits, it 
turned out that typically a somewhat smaller value of the expectation value for S was called 
for, say, around 1/2. Also the other VEVs would be somewhat changed in the "improved" fit 
including this "factorial correction" (so called because we have seen that it is square roots of 
factorials that come in). It is clear that the parameter £, which roughly plays the role of the 
Cabibbo angle and also explains the ratio of the mass of the first- to the second-generation 
(as being £ 2 ), will be smaller with the "factorial corrections" included in a fit because of the 
highly suppressed first-generation masses, of course, tended to have more fields to permute and 
thus a larger enhancement due the "factorials" than, say, the second-generation. Such a fit, 
then, must make £ smaller with "factorials" in order to compensate so as to keep the first- to 
second-generation mass ratio fixed. Indeed the £-VEV without "factorials" is fit to about 1/10 
whereas with "factorials" is instead rather 1/30. 

Here is another little detail: Since the field S has the expectation value near unity (in Planck 
units) having one or several S*-factors is not distinguishable from the fitting. Hence, S'-factors are 



W 



X 



Fig. 2 

One of the typical Feyn- 
man diagrams getting mass 
matrix element. The tad- 
pole lines with crosses sym- 
bolise VEV of the different 
Higgs fields. 
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not very well determined from phenomenology. Now the number of 5-factors needed to make 
a given quantum number shift can be changed by modifying the assignment of the quantum 
numbers for some of the other Higgs fields. Hence one can compensate for the changed quantum 
number by accompanying the other fields with a number of S-fields, so as to get the the same 
total transition in quantum number. 



5 Data and fitting of charged sectors 



Before the results are presented, we should review briefly the neutrino experimental data |l|, 
11, |l^, 13, 14, [l5|, |l6|, [l?], |l8|| : the best fit to the Super-Kamiokande atmospheric neutrino data 



shows near-maximal mixing ang le and Am 2 tm ~ 3.2 x 10" 3 eV 2 , the 90% C.L. range being 
(2 — 7) x 1CP 3 eV 2 . As for the solar neutrino data there are four different solutions; the MSW- 
SMA require values of the mass square difference and mixing angle that (at 99 % C.L.) lie in 
the intervals 

4.0 x 10~ 6 eV 2 < Am% < 1.0 x 10~ 5 eV 2 , 
1.3 x 10~ 3 < sin 2 26» Q < 1.0 x 10~ 2 
(3.25 x 10~ 4 < tan 2 6» Q < 2.5 x 10" 3 ) , 

whereas the MSW-LMA solution is realised in the intervals 

7.0 x 10~ 6 eV 2 < Ami ~ 2 -° x 10 ~ 4 ey2 . 
0.50 < sin 2 26» < 1.0 (0.17 < tan 2 6> < 1.0) . 

The LOW solution lies approximately in the region 

0.4 x 10" 7 eV 2 < Am% Z 1.5 x 10~ 7 eV 2 , 
0.80 < sin 2 26» < 1.0 (0.38 < tan 2 6 & < 1.0) , 

and the Vacuum oscillation (VO) in 

j 

0.67< sin 2 29 Q < 1 (0.27< tan 2 9 Q < 4) 



5 x 10~ n eV 2 < Am 2 Q ~ 5 x 10 ~ 10 ey2 



Since the parameter values resulting from the fits to the charged mass matrices for, e.g., £ 
are of the order of a factor 3 smaller than without the "factorial correction" , we have a priori 
an uncertainty of that order if we are not careful to treat the neutrino part of the calculation in 
the same way as the charged particle sector. In our previous paper we made partial "factorial 
corrections" for the neutrinos using the "old" parameters, of which £ is most important, taken 
from fits to charged masses and mixing angles without "factorial corrections" . One can immedi- 
ately foresee that making "factorial corrections" for both charged fermions and neutrinos should 
have the same effect as correcting the previous work j|] - very roughly - as if we simply decrease 
the £-value used for the neutrino predictions by about a factor 3. 

In several fits to the masses and mixing angles for the charged lepton and quark were 
found using different quantum number assignments for the Higgs fields^ £, T, W, <p ws with 

5 In the paper ^ the gauge group does not include the see-saw sector; therefore in fitting the neutrino masses and 
mixing angles using right-handed neutrinos, we have to add new gauge fields and Higgs fields which spontaneously 
break the U(1) B -l gauge group; i.e. we have to introduce the additional parameters, namely e and r\. They can 
also be taken in the range |e| < 1 and \r\\ < 1. See section 3 for details. 
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Table 2: Fits including 0(1) factors to quarks and charged leptons . The VEV of Higgs fields 
are measured in the Planck unit. For the notation of a , /3 ,7 , 5, see formulas (Q). 





a 


p 


7 


5 


(W) 


(T) 


(S) 


(0 


~2 

X 


i 


-1 


-1 


-1 


1 


0.0741 


0.0635 


0.487 


0.0331 


1.57 


ii 


-1 





-1 


1 


0.0945 


0.0522 


0.347 


0.0331 


1.41 


hi 


-1 





1 


-1 


0.0857 


0.0522 


0.686 


0.0365 


1.59 


iv 


-1 


1 


-1 


1 


0.0894 


0.0525 


0.756 


0.0247 


1.26 


V 


-1 


1 


1 


1 


0.0945 


0.0474 


0.653 


0.0365 


1.46 


vi 








-1 





0.0741 


0.0810 


0.286 


0.0300 


1.27 


vii 








-1 


1 


0.0857 


0.0548 


0.442 


0.0347 


1.40 


viii 














0.0816 


0.0735 


0.299 


0.0331 


1.37 


ix 





1 


1 





0.0945 


0.0522 


0.721 


0.0331 


1.62 


X 


1 


-1 


-1 


1 


0.0857 


0.0522 


0.622 


0.0300 


1.44 


xi 


1 


-1 


1 


1 


0.0900 


0.0497 


0.622 


0.0383 


1.70 


xii 


1 





-1 





0.0900 


0.0522 


0.537 


0.0422 


1.79 


xiii 


1 





-1 


1 


0.0816 


0.0549 


0.538 


0.0331 


1.64 


xiv 


1 


1 


-1 





0.1042 


0.0522 


0.346 


0.0444 


1.85 



respect to changing them by the quantum number combination of S, where we also vary (S) 
rather than fixing (S) = 1. The best fitting quantum number assignments from formula within 
the range \a\ < 1, \(3\ < 1, I7I < 1 and |<5| < 1 are reviewed in Table The notation x 2 m the 
last column of Table Q is misleading with respect to its normalisation in as for as it is defined 
as what \ 2 would be if the uncertainty in the logarithms were unity: 



X' 



E 



1 2 



in 



in 



m exp J 



(13) 



where m are the fitted charged lepton and quark masses and mixing angles and m e xp are the 
corresponding experimental values. The Yukawa matrices are calculated at the fundamental scale 
which we take to be the Planck scale. We use the first order renormalization group equations 
(RGEs) for the Standard Model to calculate the matrices at lower scales. Running masses 
are calculated in terms of the Yukawa couplings at 1 GeV. A typical result of a fit including 
averaging over the order of unity {i.e. 0(1)) random numbers (see last of the paragraph of 
section 2) is represented in Tabic |3|.p] 



6 Results and discussion 



In Table 4 we have, in addition to presenting the measured numbers for the ratio of the 
mass square differences and the mixing angles for neutrinos, given a series of predictions for 
these quantities for various quantum number assignments of the Higgs fields in our model. The 

6 Really Table 3 which is copied from Ref. is said to have 0(1) factors in distinction to another table in 
that article in which only random phases were used. 
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Table 3: Typical fit including averaging over 0(1) factors with a = — 1, /3 = 1, 7 = 1 and 5 = 1. 
All quark masses are running masses at 1 GeV except the top quark mass which is the pole 
mass. 





r iiueti 


Hz -X.pt!! llIlClllcLL 


m u 


3.1 MeV 


4 MeV 


m d 


6.6 MeV 


9 MeV 


m e 


0.76 MeV 


0.5 MeV 


m c 


1 9Q PoV 

i.zy Vjev 


1 /I PoV 

±.4 v^e v 


Ttl „ 


390 MeV 


200 MeV 


m M 


85 MeV 


105 MeV 


M t 


179 GeV 


180 GeV 


m b 


7.8 GeV 


6.3 GeV 


m T 


1.29 GeV 


1.78 GeV 




0.21 


0.22 


v cb 


0.023 


0.041 


v ub 


0.0050 


0.0035 


Jcp 


1.04 x 10" 5 


2-3.5 x 10" 5 


x 1 


1.46 





first column in this Table - below the experimental summary part - specify using a capital 
Roman number I, II or III the quantum number assignment for the B — L charge breaking field 
(j) B -L- The dominant element of the right-handed Majorana mass matrix, is for the well-fitting 
scenarios, always one from the third column (the v Rt column). The capital Roman number 
denotes then the row number of this dominant element. Which element dominates depends of 
course on the <\> B -l quantum numbers (see equations @, @ and ([l0[)). Furthermore, there is 
a small Roman number referring to the quantum number combination chosen from Table |2| for 
the Higgs fields £, T, W and <f) ws which were already used in the charged sector fitting. We only 
use combinations with good fitting of the charged sector. We also have the freedom to choose 
the quantum numbers of the fields only of relevance for the neutrinos x an d <Pb-l by shifting 
the Higgs fields S as parameterised by e and rj. We have marked the best fit - case III — ix with 
the parameter choice 77 = and e = —1 - with the black bullet in Table 4. 

The VEV of 4> B -l is irrelevant when we only look for the ratio Atuq/ Am^ tm but not for 
the absolute values of these mass square differences. Therefore we do not have to fit the VEV 
of 4> B -l- But the x-field VEV is fit - by hand - for each of the quantum number combinations 
considered and this VEV relative to the fundamental units (i.e. the suppression factor really) is 
presented in the last column. 

In spite of the fact that our previous article used only the case where (Pb—l had the quantum 
numbers of the combination u Re v RT which we here call case I, we do not, upon using the factorial 
correction fitted expectation value of £, T, W, and cj) ws (especially £) get as good fits to the 
neutrino mass matrices as before. In fact, the case I calculations, using various well-fitting 
assignments in the charged sector, consistently, yield predictions for the solar to atmospheric 
mass square difference ratio, Amg/Am^ that are too small by about one order of magnitude 
compared to the MSW-SMA phenomenological value. This is due to that we - ignoring factorial 
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Table 4: The numerical results of the three different cases, I, II and III combined with various 
charged sector fits denoted by small Roman numbers. The best fit, case III — ix with r\ = and 
e = — 1, is marked with black bullet. See more detail in the text. 





tan 2 6> 

[sin 2 26» Q ] 


A ™2 




tan 2 6> atm 

[sin 2 20atm] 


tan 2 (9 e3 

[sin 2 20 e3 ] 


(x) 


experi- 
mental 
data 


VO 


0.23 - 1.0 

[0.6-1.0] 


« 10" 7 


0.485 - 1.0 
[0.88 - 1.0] 


<0.026 
[£0.1] 


- 


SMA 


(0.33 - 2.5) x 10~ 3 
[(1.3-10)xl0" 3 ] 


1.5 to 5 7 x 10~ 3 


— 


LMA 


0.17-1.0 

[a c 1 nl 

[U.o— i.uj 


9.4 x 10~ 3 


- 


LOW 


0.38 — 1.0 

[0.8-1.0] 


O 1 +11 ^ , in — 5 

3.1 _ 2 3 x 10 




I-xiii rj = 0, e = 1 


7.8 x 10~ 4 


6.7 x 10~ 5 


1.00 


8.0 x 10" 4 


0.018 


I-xiii rj = —1, e = 1 


7.8 x 10" 4 


7.8 x 10 _b 


1.00 


7.8 x 10~ 4 


0.038 


Il-ii 77 = 1, e = 


3.8 x 10~ 3 


1.0 x 10~ 2 


0.93 


8.8 x 10" 4 


0.016 


Il-vi 77 = 1, e = 


4.7 x 10" 3 


5.8 x 10" 3 


0.96 


1.0 x 10~ 3 


0.018 


Il-ix 77 = 0, e = 1 


1.1 x 10~ 3 


1.0 x 10~ 2 


0.96 


2.1 x 10" 4 


0.026 


II-v 77 = — 1, e = 1 


1.8 x 10~ 3 


8.9 x 10~ 3 


1.00 


2.8 x 10~ 4 


0.022 


Il-ix 77 = —1, e = 1 


1.2 x 10~ 3 


9.0 x 10~ 3 


1.00 


1.9 x 10" 4 


0.022 


Il-xi 7/- = —1, e = 1 


2.1 x 10~ 3 


8.4 x 10~ 3 


0.95 


2.6 x 10" 4 


0.023 


Il-vi 77 = 1, e = — 1 


4.3 x 10~ 3 


6.6 x 10~ 3 


0.92 


1.0 x 10~ 3 


0.018 


lll-v 77 = e = 


1.1 x 10~ 3 


8.8 x 10~ 3 


0.93 


6.6 x 10~ 4 


0.032 


III-viii 77 = 0, e = 


2.9 x 10" 3 


8.7 x 10" 3 


0.99 


2.7 x 10" 4 


0.063 


Ill-xi 77 = 0, e = 


1.4 x 10" 3 


8.0 x 10" 3 


0.99 


7.3 x 10" 4 


0.030 


TTT-v 7? — f) f — 1 


1.2 x 10~ 3 


6.8 x 10~ 3 


1.00 


6.4 x 10" 4 


0.031 


III-viii 77 = 0, e = 1 


3.0 x 10~ 3 


7.8 x 10~ 3 


0.95 


2.3 x 10" 4 


0.065 


Ill-xi 77 = 0, e = 1 


1.5 x 10" 3 


6.3 x 10" 3 


0.98 


6.4 x 10" 4 


0.030 


Ill-ii rj = 1, e = 1 


3.6 x 10~ 3 


1.0 x 10~ 2 


0.92 


1.7 x 10~ 3 


0.018 


Ill-vi 77 = 1, e = 1 


4.4 x 10" 3 


8.0 x 10" 3 


0.96 


2.3 x 10~ 3 


0.019 


i T T_Tr 1 1 — 1 a — 1 

±±± V // X, c J- 


1.5 x 10~ 3 


6.3 x 10~ 3 


0.97 


6.2 x 10" 4 


0.028 


Ill-ix 77 = —1, e = 1 


1.0 x 10" 3 


6.2 x 10~ 3 


0.95 


4.2 x 10" 4 


0.028 


Ill-ii 7/ = 1, e = 


3.6 x 10~ 3 


1.0 x 10~ 2 


0.92 


2.4 x 10~ 3 


0.019 


III-v 77 = -1, e = 


1.5 x 10~ 3 


6.8 x 10~ 3 


0.94 


6.2 x 10~ 4 


0.028 


III-viii 77 = -1, e = 


3.0 x 10~ 3 


7.5 x 10~ 3 


1.00 


2.5 x 10" 4 


0.130 


Ill-ix 77 = — 1, e = 


1.0 x 10~ 3 


6.7 x 10~ 3 


0.94 


4.2 x 10" 4 


0.028 


Ill-xi 77 = —1, e = 


1.8 x 10" 3 


6.4 x 10~ 3 


0.99 


6.3 x 10" 4 


0.028 


III-v 77 = 0, e = -1 


1.2 x 10~ 3 


7.5 x 10~ 3 


0.98 


6.4 x 10" 4 


0.031 


• Ill-ix 77 = 0, e = — 1 


8.3 x 10~ 4 


5.8 x 10" 3 


0.97 


4.3 x 10" 4 


0.035 


Ill-xi 77 = 0, e = -1 


1.5 x 10" 3 


7.8 x 10~ 3 


0.97 


6.4 x 10 -4 


0.030 


III-iii 77 = 1, e = — 1 


1.5 x 10~ 3 


7.2 x 10~ 3 


0.97 


4.5 x 10~ 4 


0.025 


Ill-ix 77 = 1, e = —1 


1.3 x 10" 3 


7.5 x 10" 3 


0.97 


2.8 x 10~ 4 


0.039 


III-v 77 = —1, e = —1 


1.5 x 10~ 3 


8.2 x 10~ 3 


0.99 


6.5 x 10" 4 


0.027 


Ill-ix 77 = — 1, e = — 1 


9.8 x 10~ 4 


7.9 x 10~ 3 


1.00 


4.5 x 10" 4 


0.027 


Ill-xi 77 = — 1, e = — 1 


1.8 x 10~ 3 


7.9 x 10~ 3 


1.00 


6.4 x 10" 4 


0.027 
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0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 0.2 

(x) 

Figure 3: The numerical results in the best case, Ill-ix with r\ = and e = —1, of the ratio of 
the solar neutrino mass square difference to that for the atmospheric neutrino oscillation (thin 
solid line), and the squared tangent of the solar neutrino mixing angle (cross-dotted line), the 
atmospheric neutrino mixing angle (thick solid line) and the mixing angle 9 e 3 (dotted line) as 
a function of the VEV of x Higgs field. The results are obtained with 100, 000 order of unity 
random number selections. 



corrections and order of unity details - get £ for this ratio in the case I quantum number 
assignments. With the slightly "unfair" treatment of our previous article (using £ etc. values 
deduced without use of factorials in the charged sectors - in spite of having "factorials" in the 
neutrino part of the calculation) - we managed to get an increase by a factor 6 for this ratio, 
but that is not sufficient with the correctly fitted £-VEV. We have therefore only presented a 
few examples with case I in Table 4. 

Both case II and case III have many well-fitting quantum number assignments associated 
with them as is seen from Table 4. Generally, the tendency for both case II and case III is 
to predict the mass square difference ratio Amg/Am^ a bit too high by half an order of 
magnitude. Also the solar neutrino angle is predicted a bit too high. 

The case III fits a little better than case II because case III predicts the solar mixing angle 
a little better, i.e. to be a bit smaller]]. In all cases our predictions of the solar mixing angle is 
in the region of 1.5 x 10" 3 . It fits very well for the MSW-SMA. It is, however, not compatible 
with the other solution regions VO, LMA nor LOW. So it is necessary - unless our model is 
drastically changed - that the MSW-SMA is upheld. Otherwise our model disagrees significantly. 
Actually, this decrease in Rvalue insures that the first generation matrix elements proportional 
to £ or £ 2 have many fields and therefore get large "factorials" that must be compensated for. 
In Table 2 we show some numerical results of the most suggestive calculations: the combination 
(a = —1, j3 = 1, 7 = —1, 5 = 1), i.e., case iv, giving the very best fit to the quarks and 

7 Equally good fits to the mass square difference ratio, Am,Q / Amatm — 5.8 x 10 -3 is obtained for solutions 11 — vi 
and III — ix, but III — ix has the solar mixing angle closer to the global experimental fit Jis|], tan 2 8q ~ 5.8 x 1CP 4 . 



13 



charged leptons supplemented with the trivial case (i.e.,?] = 0, e = 0) gave no modification of 
the already rather simple quantum number combinations for x an d ^b-l in Table 1. This fitting 
of only x to neutrinos (case III — iv with r/ = and e = 0) gives tan 2 9q = 7.5 x 10~ 3 and 
Amg/Amj tm = 1.3 x 10 -2 . We also considered some variations in r\ and e in order to study 
other than the very best fits of the charged sector. The main point is that the different choices 
do not yield very different fits to the neutrino oscillations and that there are very many good 
fits of which a selection - the best mainly - has been put in the table. 

In the examples with case II and also case III, we tend to get close to the MSW-SMA 
solution although the solar neutrino mass square difference Ar/ig is predicted a bit to the high 
side of experiments, i.e. AttIq/ Am 2 tia ~ 6x 10~ 3 compared to the experimental Arrig / Am 2 tm ~ 
2 x 1CP 3 . Using the best fitting choice of which S-quantum numbers are added to the Higgs 
fields [a. = 0, j3 = 1, 7 = 1, 8 = 0, rj = 0, e = —1) gave results just of that character (see also 
Fig. 3): 

Am% +30 o . 1 

- 5.8 x 1(T 3 (14) 



+21 

tan 2 = 8.3 x 1(T 4 . (15) 
—6 

The deviation from the experimental data of our predictions should be considered within un- 
certainly limits. We interpret our "being uncertain by order unity" as meaning a relative un- 
certainty ±64% for masses and mixing angles. Note that for the mass square difference ratio 
AmQ/Am 2 tm , which involves four mass factors, the uncertainly is larger by a factor 2v2 ~ 3 
counted in logarithmic uncertainty meaning an uncertainty of the order "J^ 1 %. But for the 
squared mixing angles the uncertainty is "only" ± 2 x 64% meaning +2 j® %. With these uncer- 
tainties the predictions are in good agreement with experiment. 

6.1 Crude calculation of the mass square difference 

The results of the see-saw mechanism is that the light neutrino masses are quadratic in the 
Dirac masses and inversely proportional to the heavy right-handed Majorana masses: 

M eS ^M^M^(M^) T . (16) 

The mass square difference ratio, Am 2 ^/ Am 2 tm , from equations ([?]), fli~0| ) and (|l~6|) can be 
approximately calculated (estimating by rules like A a a + Aft b ~ \J a 2 + b 2 , where A a and Aft are 
order of unity random numbers) in the best case, i.e., Ill — ix with rj = and e = — 1, as follows: 

part of 



U \2£ 2W 2 T22( , )2 /W35ST£ 2 6(K/l4S 3 r£ 3 
M eff , casc „, ix ~ (^s) b WTx (gg-j) g2 Ti 2y /3 T 

l " = °" = - 1 detM * W70S 2 £x 2^6 X 



33/2 2V6£ 3\/21/2xA / 6^355T^ 2 6\Z35S 2 Tf 
2^6^ V^46£, 2 I2y/$xi 2 60VUS 3 T£ 3 2y/?,T 



^3^21/2 xi 12^5 xi 2 3^/165 x 2 £, 2 J \60Vl54 S 3 T f X ^VTETx 




part of M^ 1 part of (MP) 1 
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<<w 2 s 2 w 2 r 2 x 2 e 



detM R 



,0 ~0 ~0 \ 

-0 1103 T 2 1391 T X 
\~0 1391 Tx 2206 x 2 J 



(17) 



where we have used the Higgs fields S, W, T, £, X and the fields S*, W\ T\ ft, x ] (with 
opposite quantum charges) equivalently because the non-supersymmetric model is considered 
here. Moreover we have used the numerical value S 2 ~ 1/2 and the relation T ~ % to obtain 
the large atmospheric mixing angle. It is important that the dominant matrix elements are not 
uncorrelated. The determinant of the dominant 2-by-2 subgroups of M c g- must be calculated 
by multiplying the 2x2 determinant of the relevant 2-by-2 sub-matrices of the three matrices, 
M?, Mft 1 and (M?) T : 



detM eff 



case Ill-ix 
T) = 0, e = -1 
(2, 3) X (2, 3) block 



(2,3) X (1,2) block detM R I (i, 



(TVM cfr ) : 



case Ill-ix 
T[ = 0, e = -1 
(2, 3) X (2, 3) block 



iJl I (1, 2) x (1, 2) block 
li) = 0,e = -l l?7 = 0,e = -l 



x det((M i f ) T ), (1 , 2) x (2 , 3 ) wock 

I 7) = 0, e = -1 

3X10 (d^Mn)' 2 

fi v in e g4^4 T 4 x 8 e 4 (Vs) 4 {(j)B _ L y 

{d^Mnf 2 



(18) 



(19) 
(20) 



So we get 



' Am| 



detM cfr 



^4 

T 



(TrM cfr ) 2 2 " 8 (21) 

Actually, the terms coming from the (1, l)-component of the M^ 1 are somewhat dominant due 
to the factorial corrections, so that the dominant 2-by-2 subgroup, the (2,3) x (2,3) of M c g, is 
to first approximation degenerate. 



7 Conclusion and resume 



We have considered a series of extensions so as to include of neutrinos of the earlier Anti- 
GUT fit to the charged quarks and lepton masses by introduction of two more Higgs field, x 
and 4>b-l- 

To fit the atmospheric neutrino oscillation experimental data it is necessary to fix the VEV 
of the Higgs field, (x) , to be in first approximation X ~ T in order to arrange the atmospheric 
mixing angle to be of order unity. We should therefore not consider the agreement of the 
atmospheric mixing angle 9 a tm as one of the significant predictions of this model: we have, 
however, managed, at least, by such a fitting with %, to obtain a large mixing angle; 9 a t m , i.e., 
the atmospheric mixing angle is essentially an input parameter in this model. 

The success of this model should be rather judged from the following predictions: 

1 ) the solar mixing angle which comes out proportional to £ which in turn deviates only by 
"factorials" and order of unity from the Cabibbo angle. 
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2) the mass square difference ratio Arag/Am^ tm . 

3) the electron- and tau-neutrino mixing angle being small enough not to be in conflict with 
the CHOOZ measurements (U e z = sin# e 3 < 0.16). 

The first of these predictions is the major reason that our model in the present form (i.e. 
with the present quantum number assignments only varied within the limits considered in this 
article) is only compatible with the MSW-SMA solar neutrino solution. The smallness of the 
Cabibbo angle results in our solar mixing angle being so small that it would stress the model 
drastically to seek to fit one of the series of large solar mixing angle fitting regions such as LMA, 
LOW or VO. 

Within the theoretical uncertainly inherent in our predictions being only of order unity - 
even if taken as specified in Ref. Q with 64% being one standard deviation - our model fits 
the data perfectly provided that we allow the MSW-SMA solution as the (possibly) valid one 
in spite of the fact that the day-night effect disfavour this solution by two standard deviation: 
Extracting from the best fit, III — ix with r] = , e = —1, the result may be stated as 

tanfl « e~ o#c , (22) 



Am atm 13 



(23) 



The collective fit of both the neutrinos and the charged sector with 6 fitted VEVs, one of 
which (S) is still close to unity and with the (4> W s) determined from the weak interaction Fermi 
constant, is so good order of magnitude-wise that it should perhaps be considered perfect. 

Since our mass eigenvalues are not tightly degenerate, it is not likely - or rather it is impos- 
sible - that the renormalization group running should be very important as corrections to our 
predictions. In fact we expect them to have almost no influence to our order of unity accuracy 
for mass ratios of particles within the same group such as, say, the (left) neutrinos. In fact, in 
the charged sector, the renormalization group mainly just corrects the masses of the quarks by 
a factor 3 or so compared to the corresponding charged leptons. 

One could perhaps complain though that we have had too many relatively complicated quan- 
tum number assignments which, although only discrete choices, may be too many possibilities 
to make our fit convincing. The Anti-GUT model with its very many gauge fields can also be 
complained about as being too complicated, but here it should be pointed out that indeed it 
should rather gain its reason for being considered by being the largest gauge group transforming 
the Standard Model fermions plus the see-saw neutrinos non-trivially among themselves without 
unifying irreducible representations of the Standard Model. 



7.1 Quantum number system of fitting results 



Presumably the honest motivation for a fit as in the present work is that we have a major 
part of the fitting going on by fitting the discrete representation choice for the Higgs fields. 
The best fit has its Higgs field (abelian) quantum numbers listed in Table 5. This Table, in 
this sense, represents quantum numbers that are derived/inspired from experiment. One can 
then imagine using this inspiration to look for regularities suggestive of the model beyond the 
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Table 5: All U(l) quantum charges of the best fit in extended Anti-GUT model. The VEV are 
presented in the Planck unit. Non-abelian representations are given by a rule from the abelian 
ones, see after Eq. (2). 



+>ws 
S 
W 

T 

X 



SMGi 



SMG 2 





l 

6 



SMG 3 






2 
3 



B-L,2 



u B _, 




1 

3 




1 

2 



VEV 



~ nr 1 * 

0.721 
0.0945 
0.0331 
0.0522 
0.0345 
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Standard Model. It is easy to estimate the VEV for 4> B -l needed in the fit by using, say, the 
example of section 6.1. Inserting 

detM R ^72S 3 x 4 f(^B- L ) 3 (24) 

into equation (19) yields 

AmL - (TrM cff ) 2 * 1.2 x 10 3 ^J^f • (25) 

^ Wb-l) 

Then from equation (|24|) together with the atmospheric neutrino oscillation data we get 

{4>b-l) ~ 6.7 x 10 11 GeV . (26) 

We present here the order of magnitude of the right-handed Major ana neutrino masses from 
the mass matrix, Mr, of the best fitting case, III — ix with 77 = and e = — 1 including factorial 
corrections: 

M Rui rs 3 x 10 6 GeV , 

Mr„ « 4 x 10 8 GeV , (27) 



Mr » 1.4 x 10 11 GeV . 



7.2 Theoretical lesson 



It is highly remarkable that we obtain the rather good fit of the mass square difference ratio, 
Am^/Am^tm; i n the two well- fitting cases II and III: were it not for the, in principle, order unity 
factors coming out of this ratio, it should have simply been unity: Amg/Ara^ m ~ 1 since the 
Higgs VEV's drop out except for the Higgs field S. The remaining in principle order unity ratio 
comes out of the way the various numbers are combined in the calculations of inverse matrices 
and matrix products and the "factorial corrections" which are due to diagram counting and the 
field S. 

All of the numbers that are not of order unity in our model, namely the Higgs field VEV (in 
Planck units) £, % . . ., except S which has anyway (5) as 1, drop out of this ratio. This is the 
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way that our model in cases II and III solves one of the at first sight almost hopeless problems 
with data: how, in fitting the large atmospheric mixing angle, to get the hierarchically looking 
large mass square difference ratio between the middle and the heaviest eigenmass neutrinos when 
the highest mass has to mix strongly (of order unity) with other eigenstates. 

It is, however, indeed possible to interpret this hierarchical mass ratio as really due to a 
hierarchical small VEV as case I does, where this ratio goes as £ 4 . However, case I, does not fit 
well precisely because it makes this ratio too small - too hierarchical. 



7.3 Our qualitative predictions for future neutrino physics 



We summarise here our predictions for future neutrino experiments assuming that this model 
were/is right: 

1) It should be SMA-MSW that becomes the best fit (the day-night effect should be tested 
by the SNO detector JH 



2) It is good that the atmospheric mixing angle is of order unity but it should not turn out 
to be very close to tan 2 atm « 1, i.e. tan 2 6> a t m = 1 + O (lower order). There are namely 
models of a similar spirit (but not our most promising fits) which behave this way - for 
example when a pair of off-diagonal elements of M e g dominate. 

3) KamLAND should not be able to see any neutrino oscillation because of SMA-MSW. 

4) Since this model only contains three light neutrinos, there is no place for the LSND ef- 



fect n 



It is important to confront our model with baryon number production in Big Bang which in 
models like ours with pure Standard Model at the weak scale must get the baryon number as 
a B — L contribution from some other scale. This model is of the see-saw type with the scale 
given by the neutrino oscillation scale. After crude estimation of the Baryogenesis, we found 
that this model predicts the right range as Ref. |^T|, |2[ studied. Since we do not have SUSY 
in our model - preferably at least - we do not have to worry about gravitino problems ifffi] . We 
can say that they simply do not exist at all or only at such high energy scales that they are 
totally irrelevant. Another potential problem |24|, is the wash out of the B — L excess. It 
is diminished compared to the typical lightest see-saw masses considered by [E^] because of our 
relatively low lightest see-saw mass, 3 x 10 6 GeV (see equation (E~ 



7.4 The "efficiency" of our model 



The number of measured quantities, which are predictable with this model, quark and charged 
lepton masses and their mixing angles (three Cabibbo-Kobayashi-Maskawa mixing matrix an- 
gles) containing Jarlskog triangle area, J C p and also two mass square differences and the three 
mixing angles for the neutrino oscillation, Baryogenesis and neutrinoless double beta decay pq] , 
is 20. 

Our model predicted successfully all these quantities using only six parameters^ which we 
8 The VEV of Weinberg-Salam we have not counted as a parameter because of its relation to the Fermi constant. 



18 



fit - VEVs of additional Higgs fields: the genuine number of the predicted parameters is thus 
14. However we have taken into the predictions two quantities namely, tan 2 9 e 3 and "effective" 
Majorana neutrino mass for which only experimental upper bounds exist. However we must 
emphasise again that strictly speaking we have almost "infinite" input parameters but since we 
only care for order of magnitude results, we should not count couplings being order of unity. 

The assumptions of quantum numbers are a bit arbitrary or fitted rather and a bit com- 
plicated although still discrete. But otherwise it is a very reasonable and expected type of 
assumptions that are used, e.g., the major assumption of all couplings at "Planck scale" being 
of order unity is philosophically one of the easiest to get. It is (almost) just use of usual dimen- 
sionality argumentation. That the only deviation from everything being of order unity comes 
from the VEVs may find a bit of support in the well-known fact that in super-conductors, VEVs 
tend often to be very small on the atomic physics scale expected a priori. 
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